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CHAPTER ONE 


INTRODUCTION 

-j 

In this era of state- to- state chemistry, it is possible 
to study elementary chemical reaction rates from reagents in 
selected initial states to products in speci fic final states, 
thanks to the molecular beam, laser excitation, laser^induced 
fluorescence, chemiluminescence and related techniques. It is 

also possible to select the relative orientation of the 

2 3 

reactants and to monitor the product alignment. Therefore,. 

it becomes necessary to use detailed dynamical theories to 

4 ' 

predict the observables. Although quantum mechanics is needed 
to describe the molecular collision events* exact ,(i . e. , converg- 
ed) quantal calculations are extremely difficult and time- 
consuming and therefore have been reported only for tuuo 
reactions*^ 

H + Hg — ■ — > H 2 + H (R1) 

and f + H 2 FH> H (R2) 

. .0 

Quasi classical trajectory (QCT) technique, * on the other hand 



2 


is easier and has been used extensively for studying the 

dynamics of a variety of model and real systems in the last 

tu)o decades. By comparison udth the available quantal and 

4 

experimental results, it has been established that, on the 

average the QCT method is a reliable tool provided purely 

quantal effects like tunneling, resonance etc. are not import 

tant. It has been particularly valuable in correlating the 

features of the potential- energy surface (PES) and the colli- 
9 

sional outcome. For example, it has been shouin that for 
substantially endothermic reactions (AH° ^ 10 kcal mol” ) of 
the type 

A + BC ^ AB + C , (R3) 

the saddle point is in the exit channel and that the reagent 

vibrational energy (V/) is much more efficient than relative 

i Q 

translational energy (T) in causing the reaction, in quali- 

1112 

tative agreement ujith the availablB experimental results. * 

1 3 

Polanyi and Sathyamurthy have refined this idea further. The 
effect of V on the reaction cross section (Sj.) is much more on 
a sudden surface (where the change in the potential from the 
reactants to products is sudden) than on a gradual surface 
(on which the change is gradual). 

In this thesis we report the results of our QCT investi- 
gation of the dynamics of the reaction, 

^ LiF + H (R4) 


Li + FH 
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luhich serves as a prototype for the family of alkali- hydrogen 
halide reactions and also a model six- center exchange reaction 

1 4 

Since the first successful molecular beam investigation 
of the reaction 

K + BrH ^ KBr + H (R6) 

the alkali- hydrogen halide reactions have received considerable 

15 

attention from several research groups. Brooks and coujorkers 
showed that for the reaction 

K + CIH — ^ KCl + H (R7) 

HCl (v = 1) is nearly 100 times more effective than ,HCl (v =*0) 

and that \I is much more efficient than T in causing the reaction* 

16 

Heismann and Loesch have reported recently that for the 
reaction 

K -f XH ^ KX + H (X =C1, F) (RB) 

U is more efficient in promoting the more endothermic reaction 
(XsF) and that the efficiency decreases at higher collision 
energies. The chemiluminescence depletion study of the sodium 
analog 

» NaX + H <X » C1, .F) (R9) 


Na + XH 
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has shouin that u/ith ouer 90/jf of the total energy (E) present 

as vibration, the collision efficiency is approximately unity. 

Similar vibrational enhancement has been observed for the 

18 

alkaline- earth analogues of these reactions also. On the 
basis of these results* these reactions are expected to have a 
sudden typ e P ES . 

When compared to the effect of T and U, the role of rota- 
tional energy (R) on reaction dynamics is not well understood. 

The effect of R on various reaction attributes has been 

1 9 

reviewed recently. There are some experimental results availa- 
ble on the effect of R on for a few alkali-hydrogen halide 

reactions. From a chemiluminescent depletion study of the 

2 0 

reaction (R9), Blackwell et al, found that there was an initial 
decrease followed by an increase in the reaction rate coeffi- 
cient (k) with increase in R, with a minimum occurring at a 
rotational quantum number, 3 1 -10. An initial decline in 

Sj,(3) has been reported * for the potassium analog and there 
are some indications of a subsequent increase in 3^(3), 

22 

Lee and coworkers have recently studied the reaction 
(R4) using the crossed molecular beam technique. From the near 
forward- backward symmetry of the product angular distribution 
(pad) at T s= 3*0 kcal mol , they concluded that this reaction 
proceeded via complex formation, with coplanar geometries being 

-1 ’ 

preferred. However, at T =8,7 kcal mol the PAD became more forward 
peaked for Li^ indicating the diminished importance of tho complex 
formation at higher collision onergios. Bocauso of tho possibili 
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of tunnelling, they were not able to ascertain the classical 
barrier height for the reaction. The effect of R and U 

on Sj, remain to be investigated for this reaction. 

Despite the number of experimental studies for a variety 

of alkali- hydrogen halide reactions the results have been inter- 

2' 3 ' ' ' ' ' 

preted mainly in terms of the simple "harpoon mechanism". 

Euan though this has been successful in explaining the large 

magnitude of and its change with, change in the metal and the 

halogen atom involved, it cannot explain certain aspects of the 

dynamics like specific energy disposal. Also, vibrational 

enhancement of S^, cannot be accounted for satisfactorily. 

24 

Polanyi and cou/orkers had proposed a simple DIPR (Direct 
Interaction with Product Repulsion) model, which provides a 
qualitative explanation for the product energy and angular dis- 
tribution, for reactions forming an ionic bond. However, this 
model is applicable for direct reactions only. 

QCT studies of the dynamics of the alkali-hydrogen halide 

exchange reactions have not been reported until recently as 

there were no accurate PES available. The only system amenable 

to elaborate ab initio investigation is the reaction (R4), The 

nonreactive part of the PES for this system was published by 

25 

Lester and Krauss which showed the presence of an attractive 
well. This has been confirmed subsequently by Trenary and 
Schaefer. Balint-Kurti and Yardley (BKY) published an 
ab initio PES including a send, empirical correction for the 
reaction in the collinear LiFH configuration. They also found 
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an attractive uiell in the entrance channel and calculated Ej^ 

- 1 

to be 21.6 kcal mol . In addition, they found E|^ to decrease 
as the LiFH angle decreased from 180°, Zeiri and Shapiro (ZS) 
reported^® a semi- empirical PES for a series of alkali- hydrogen 
halide reactions including (R4) , They predicted the lomest ualue 

of Ej^ to be 12.4 kcal mol”'^ at 0(LiFH) = 105°. An extensive 

ab initio PES for (R4) was published recently by Chen and 
29 

Schaefer. Their PES for the collinaar geometry mas found to 
be similar to that of BKY in many respects. The classical 
endoergicity was calculated to be 2.9 kcal mol , but with 

zero-point corrections, the reaction mas found to be exothermic 

- 1 

by 1.7 kcal mol . Their surface also contained a mail in the 

- 1 

entrance channel mith a mell depth of 4,5 kcal mol corres- 
ponding to the bent complex mith S( LiFH) =1 14° , The saddle point 

- 1 

mas located in the exit channel, 10 kcal mol above the reac- 
tants mith S( LiFH) =7 4°. This jsur face has been fitted to an 

30 

analytic function by Carter and Murrell, after appropriate 
scaling as suggested by Chen and Schaefer. The standard devia- 
tion of this fitted surface (CSCM surface) from that of the 

. m.'j 

scaled ab initio points mas 1.9 kcal mol mith a maximum error 

- 1 

of 6 .7 kcal mol . 

After this thesis mork mas begun, some quantum mechanical 

(QM) as mell as QCT results on ZS surface mere reported in the 

literature. In their QM study of the collinear reaction on the 

31' 

ZS surface, Walker at al. found a large number of resonance 

. " R 

peaks in the reaction probability (P ) E curves suggesting 



7 


the formation of long-lived complexes. For T>11.5 kcal mol 

the HF molecule from its ground vibrational state mas found 

R 

to react mith significant P , , Homever, the collinear QM calcu- 
32 

lation on the CSCM surface resulted in a negligible value of 

0.05) for HF (v = O) at T = 21 kcal mol“^. On the other 

R 

hand, mhen \i mas higher than P mas nearly unity at very 

33 

lorn values of T. Zeiri et al. have reported the results of 
their QCT calculations on their (ZS) semi- empirical PES. Their 
calculated S of 2,46 % and the pronounced backmard scatter- 
ing of the product molecule mere not in agreement mith the 
experimental result meaning that thi^ surface is not accurate 
enough for describing the LiFH dynamics. Therefore it becomes 
necessary to investigate the dynamics of this prototype alkali- 
hydrogen halide exchange reaction on an accurate surface. 
Fortunately, a reliable ab initio PES (CSCM) is available for 
the collinear as mell as noncollinear geometries and me have 
carried out detailed QCT studies for this reaction on this 
surface. We compare our results mith the available experimental 
results and assess the accuracy of the ^ initio PES. A detailed 
investigation of the effect of V and R (in comparison to that 
of T) on for this reaction mould shed light on the dynamical 
factors operative in this as mell as other members of this 
family of reactions. 

We report in Chapter 2 the results of collinear QCT 
studies for this reaction. Details of the three dimensional(3D) 
dynamical studies and results are given in Chapter 3, mith 
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special attention being given to the effect of R and the 
relative orient ation, on various reaction attributes for this 
reaction. 

In this thesis, ive also report on the development of a 
softiuare package to study the detailed dynamics of a six-centre 
(6-C) exchange reaction using the QCT technique and present the 
preliminary results obtained for a model system. 

Historicslly, the reaction 

H 2 + I 2 2 HI (RIO) 

had been considered, to proceed in one step through a four- 

centered (4-C) transition state. In the sixties, it mas shomn 

34 35 

to proceed via a chain mechanism involving atomic iodine. * 

Raff et al.^^ had subsequently examined the dynamics of this 

37 

reaction on a semiempirical PES by the QCT technique and 

concluded that the direct reaction of H 2 with I 2 was prohibited 

by dynamical effects and that the reaction occurred primarily 

via the formation of an H 2 I complex. Their conclusions were 

corroborated by the combined phase-space/trajectory calcula- 

3 8 

tions of Daffe et al. By a molecular beam investigation 

39 

using a supersonic nozzle, Daffe and Anderson had earlier 
found the reverse reaction 

HI + DI » HD + I2 , (RTl) 

not to occur even with a relative translational energy of 
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- 1 

upto 110 kcal mol . They suspected that vibration or rotation- 
al excitation of the co-Ll-iling partners u/as required for the 
reaction to proceed through the molecular mechanism. Recent 

4q 

experiments by Horiguchi and Tsuchiya shouj the vibrational 

excitation of the HI molecules to enhance the rate of this 

reaction. The other potential 4-C reactions which have been 

414 2 

studied by the QCT technique are * 


+ 

CM 

X 

>^2 ~ 

— ^ 2 HF 

(R12) 

HBr + 


— » HCl + BrCl 

(R13) 

HBr + 

BrCl — 

— ^ HCI + Br2 

(R14) 


and the results indicate that these reactions also must be 
proceeding mainly via chain mechanism. 

The hydrog en- deuterium exchange reaction 

H 2 + D 2 ^ 2 HD (R15) 

has been in the forefront because of the severe conflict botiueen 

43 

the theoretical and experimental results/ Shook-tuha studies 

44 45 46 

of Bauer and Ossa and others * yielded an activation 
energy (E^^) of 42 kcal mol . The results were explained by 
postulating a vibrational excitation mechanism via a 4-C tran- 
sition state. Initial doubts on the reliability of the shock*- 

tube results were removed by the stimulated Raman experiments 

4V 40 49 

of Bauer and cpworkers* Yet several inv/estigators * maintain 

that only the atomic mechanism can sufficiently explain all 
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the observations. Also, a recent laser^beam investigation^^ 
of the reaction 

HD + HD > H 2 + D 2 (R16) 

shouied that this bimolecular reaction did not occur even u/han 
one of the HD molecules o/as excited upto v = 5 level. 

These intriguing experimental results had led to a series 

of theoretical investigations of the reaction (R15), An exten- 

51 

sive initio search of Silver and Stevens failed to 
produce a 4- C transition state, luith an loiuer than the disso- 
ciation energy of the hydrogen molecule. The collinear symme- 
trical transition state, ujhich lies below the dissociation 
limit for H 2 , cannot serve as the transition state for a 4-C 
exchange reaction since the end atoms must depart from the 

H^/Dj 

52 

collisions was carried out -by Brown and Silver on a London- 
type PES which showed that only-N^Tjlf of the total number of 
trajectories led to the products when both the molecules were 
in the ground vibrational state (v - O) , However, when both 
the molecules were excited to v »3 level no reaction was observed. 

53 

Another study on a semi- empirical repulsive type PES, however, 
resulted in rv 3 / reaction when both the molecules were excited 
to V s 3 . 

The failure to realise a 4-C transition state has been 

54 

explained on the basis of Woodward-Hof fmann rules. Since 


collinear array to recombine. A detai^led QCT study of 
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the (electronically) ground state reactants are not correlated 
u/ith the ground state of the products because of the conserva- 
tion of orbital symmetry, a large barrier for the reaction 
(R15) results. This led Wright^^ to suggest a 6- C transition 
state as an alternative. His LCAO-MO-SCF calculations^^ yielded 

- "I 

an Ej^ of 90 kcal mol , ujhich is ujell below the dissociation 

energy of H^, for the formation of from 3 H 2 . Using London- 

57 

Eyring-Polanyi formalism, Thompson and Suzukawa predicted a 

- 1 

value of Ej^ = 88 kcal mol whereas a large basis set ^ initio 

5 8 .pi 

calculation of Dixon et al. yielded an Ej^ of 69 kcal mol” . 

The arguments in favour of the 6-C transition state gained 

further support, through the crossed molecular beam experiments 

of Herschbach and coworkers^.^*^^ They found that the reaction 

between CI 2 and Br 2 or HI did not occur at collision energies 

- 1 

as high as 25 kcal mol . However, at high pressures and 
low temperatures in the supersonic nozzle expansion where there 
was a high concentration of ( 012)2 they found the reaction to 
proceed at collision energies as low as<v3 kcal mol” . This 
they ascribed to the 6-C reactions: 


(012)2 ^ ^ 

( 012)2 + ^ HCl + ICl + CI 2 (RIB) 

A similar 6-0 reaction 


(HI)2+ F 2 


HF + HI + IF 


(R19) 
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has been observed by Durana and McDonald. 

Therefore, it would be worth investigating the detailed 
dynamics of a 6-C exchange reaction by the QCT method on a 
realistic PES. Such a study would also enhance our understand* 
ing of other interesting processes like collision-induced disso- 
ciation of a van der Waals molecule, and vibrotational energy 
transfer of a diatom by a van der Waals species. The details 
of the development of a software package (for thb first time) 
to study a model 6-C exchange reaction along with some preli- 
minary results are presented in Chapter 4, Summary and conclu- 
sions of the thesis are given in Chapter 5, 



CHAPTER TWO 


DYNAMICS OF THE COLLINEAR Li -f FH — ^ LiF + H REACTION 

2,1 Computational Procedure 

27 

Balint-Kurti and Yardley have reported the _ab initio 

potential energy (PE) values over a rectangular grid of R,j = 

^LiF “ 2.0, 2.5, 3.0, 3.2, 3,5, 3.8, 4.0, 4.5, 5.0 and 6.0 au 

(1 au = 0.529167 %) and R 2 = ^hF 1 .25, 1 .50, 1.76, 2.00, 

2 ,25, 2.50, 2.75, 3.00, 3.50, 4.00, 4.50 and 5.00 au. Houuever, 

for QCT calculations, PE values are needed at larger R^^ , i *1,2 

values. We have assumed the interaction potential to become 

zero at Rj^ *= 10,0 au. The LiFH potential at this distance mould, 

62 

therefore, be the initio value for the isolated LiF and 

HF molecules. The PE values at R^ =7,8 and 9 au and R 2 = 6, 

7, fl and 9 au mere calculated using an R’"'^ interpolant. The PE 

and its first derivatives mith respect to R-j and R 2 et an 

arbitrary geometry mere calculated by the tmo-dimensional spline 

63 64 

interpolation procedure. * The resulting PES is displayed 

65 

in the form of a contour plot on scaled and skemed coordinates 

in Fig. 1. The ab initio (CSCM) surface of Chen and Shaefer 

30 

as fitted by Carter and Murrell is plotted in Fig. 2, The 
functional form and the parameters used are given in Chapter 3. 
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We have constructed a semi- empirical London-Eyring-Polanyi-Sato^ 
(LEPS) surface u/ith E^ equal to that of the BKY surface and it 
IS shown in Fig. 3. The functional form and the parameters for 
this surface are given in Appendix I. We have carried out 
collinear QCT calculations on all the three (BKY, CSCM and LEPS) 
sees. The details of the calculations can be found else- 
uihere.'^ For a fixed initial U and T, 50 trajectories were run, 
varying the vibrational phase systematically between 0 and 271 , 
ujith all the trajectories starting at 12 au on the BKY and 
LEPS surfaces. On the CSCM surface, we had to start the trajec- 
tories at = 22 au as the potential was not negligible at 
shorter distances, 

-Besults and Disnufisinn 
R 

The P values, calculated using the BKY and CSCM surfaces - 
are plotted as a function of \1 at different values of E = T y = 
1.5, 2.0, 2.5, and 3.0 eV , in Fig. 4. Both surfaces lead to 
nearly the same result. Even at the highest E = 3.0 eU, P*^ is 
zero for U corresponding to vibrational quantum numbers v = 0 
and 1, although T is at least two times E^. However, if v is 
changed to 2, P^ rises to unity even at low T. In order to 
assess quantitatively the vibrational threshold, we varied U 
in small intervals near \l = 1.0 eV, at E = 1.5 eU. it is found 

anaray required to cross the barrier is in 
or reagent vibration, the reaction does not take place . 
There is very little effect of T on P^ on both surfaces as 



Pl^ faujr 



RufCou) 

. . . f ■ 

Fig. 3 Same as fig.t for th« LEPS surface. 
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shown in Fig. 5 for v==2. This is confirmed by the QM results 
of Alvarino et al.^^ on the CSCM surface. The iiM results also 
show a marked vibrational enhancement noticed in the QCT 
results except for the nonzero due to tunnelling in the 
former. Except for the oscillations, there is, a good agreement 
between the two results as shown in Fig. 5. Similar agreement 
has been obtained for other systems.^® 

We have investigated the product vibrational energy dis- 
tribution on the two surfaces for v = 2 at different T. The 
QCT state-to-state transition probabilities displayed in Fig, 6 
for the GSCM surface show a bimodal (double-peaked) distribu- 
tion in v' at low T. With increase in T the distribution 
shifts to higher v* and loses the bimodal character. These 
results are in qualitative agreement with the QM results* 
Unfortunately, we could not make a quantitative comparison as 
the two results are not available at identical energies. 

Product vibrational energy distribution on the BKY surface 
does not show any significant bimodality (Fig. 7) but shifts 
to higher v’ with increase in T . The T dependence of the 
average fraction of the product vibrational energy ( <(f^:) ) is 
identical on both surfaces at low T as shown in Fig. 8. At 
higher T, however, CSCM surface leads to a larger <(fy/> • 

The vibrational threshold equal to noticed for the 
reaction (R4; should bo attributed to the ,§uddjgn nature of the 
PES ^ arising possibly from a change in the bonding character- 
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tM) 

Fig . ;S Comparitofi of on th« CSCM syrta c« * 
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VIB. QUANTUM NUMER OF LiF. 
Fig . 7 Same as fig. 6 for the BKV surface . 



Ftg . 8 Comporisofi of 
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covalent in the reactant (HF) to ionic in the product ( Li F) . 

As the reagents approach ea©h other, there is very little 

increase in the PE along the approach coordinate and then there 

is a sudden increase in the PE toiuards the barrier crest. Since 

most of the barrier is located in the product channel, there 

is a need for almost all of the energy required to cross the 

barrier to be in the form of W. The presence of a negligible 

fraction of the barrier in the entry channel accounts for the 

near-nonexistence of a translational threshold. This is 

perhaps the first time that such a threshold has been predicted. 

It would be interesting to find out whether such a threshold 

persists in 3D collisions. The chemiluminescence depletion 
. 17 

experiments of Bartoszek et al. showed that the sodium analog 
of the reaction (R4) does not occur unless \I is greater than Ej^ 
and under these conditions the collision efficiency is near 
unity, even if T is as low as 4/ of the total energy. Therefore, 
we infer that the reaction (R9) has a sudden type PES. 

So far our discussion has been restricted to the dynamics 
on an initio PES. For most of the systems, ab initio 
surfaces are not available and are not likely to be available 
for some time to come because of the complexity of the electro- 
nic structure calculations for many- electron systems. Therefore, 
the tendency has been to use semiempirical or even empirical 
surfaces with spectroscopic parameters for the reactants and 
products as the only input. LEP3 surfaces have been used in 
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the past for a variety of systems if only to identify the major 

qualitative features of the PES by predicting product energy 

distribution and comparing the results with the IR chemilumi- 
69 

nescence data. There are several limitations to this func- 

70 

tional form. Schor et al, have shown that the LEPS function 
could not reproduce the sudden nature of the PES for the reac- 
tion 

Be + FH ^ BeF + H (R20) 

This is not surprising as the LEPS formalism was developed 
originally for H + H 2 type systems. By searching a three 
dimensional grid of sato parameters in the range -14 0 <-1, we 
formulated an LEPS surface that met various criteriaS (i) saddle 
point in the exit channel, (2) Ej^ same as that of the BKY surface, 
(3) not too unrealistic potential well.or barrier in the 
entrance/exit channel and (4) the overall nature of the surface 
is sudden. In order to test its utility, we have computed the 
trajectories on this surface and the results are shown in Fig. 9 . 
The reaction takes place for all v and there is no vibrational 
threshold. There is, however, a translational threshold for 
V = 0 and 1 as shown in Fig, 10. This arises from the fact 
that the LEPS surface is not as sudden as the other two surfaces. 
We could not get an LEPS surface which had more sudden character 
without simultaneously introducing an unrealistically deep 
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potential uiell ($.10 kcal mol“^) in the entry channel. 
Therefore, tre had to abandon the idea of developing a model 
LEPS surface that would be suitable for a detailed dynamical 
investigation of the reaction (R4) and its analogs. 



CHAPTER THREE 


LiPH DYNAMICS IN THREE DIMENSIONS 

3.1 Pot snti al~ Energy Surface 

We had used three different (BKY, CSCM and LEPS) surfaces 

for studying the collinear collisions of the reaction (R4) but 

only one of them, the CSCM surface has been used for studying 

the 3D collisions. The BKY surface iras restricted to collinear 

27 

geometries. Although Balint-Kurti and Yardley published PE 
values for a f eui ngncollinear geometries and investigated the 

dependence of E^^ on the LiPH angle, there u/as not enough infor- 
mation for a numerical interpolation of the surface nor could 
any reasonable analytic function be fitted to the data to yield 
a PES for the noncollinear geometries. The LEPS surface has 
provision for the noncollinear geometries but does not lead to 
results similar to that of the two ab initio surfaces available 
for this system. Therefore, it is not expected to be realistic 
for the noncollinear geometries. Also, the LEPS formalism leads 
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generally to the collinear conformation as the most, preferred 
pathway whereas the ab initio calculations clearly, favour side- 
ways approach* As a result, we have not used the LEPS surface 
in our 3D dynamical investigation of the reaction (R4), 

The CSCM surface^^ is defined by the analytical functions 

^(^ 1 * r2f r^) = l/^(r^) + ^2(^2^ + ^3(1^) + ^ 2 * ^3^ 

uihere the subscripts 1, 2 and 3 stand for the molecules LiF, 

HF and LiH respectively. The diatomic potential is given 
by 

U^(r) =”Dg(l + a^p+ 32 P^+ aj ) exp (-a, P ) 

where P= r-fg* The triatomic potential is given by 

3 

^LiFH ^2* ^3^ [1-tanh (Yj^Pi/2)J 

X \J j ( 1 + 2C.p. + S C..p. P- + •*.*) 
i i^ j 

The values of. the parameters are listed in Table ,1. We describe 
here the different aspects of the CSCM surface and point out its 
shortcomings also. Although we had displayed the PE contours 
for the collinear configuration in scaled and skewed coordinates 
in Fig, 2, we reproduce the same in the (r^, r 2 ) coordinates 



Parameters for the CSCM surface 
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in Fig* 11 for comparison luith the PE contours for the LiFH =114° 

& 74° in Figi 12 & 13 respectively. The value of Ej^ decreases 

1 * • 1 ' ’ 

from 20.5 kcal mol to 7,8 kcal mol as the LiFH is changed 
from 180° to 74° and then it increases sharply as the angle is 
decreased further. The saddle point is located in the exit 
channel for all the three different configurations but the 
sudden nature of the surface decreases as the LiFH angle is 
decreased from 180° suggesting that the vibrational threshold 
noted for the collinear geometry (Chap, 2) is not likely to 
persist in the bent collisions. There are also other changes 
in the PES as the LiFH changes. The depth of the uie.ll in the 
entry channel changes as can be seen from Table 2. During a 
detailed investigation of the features of the PES, me have 
discovered that there are ‘unexpected’ additional hills and 
valleys in the pot enti al- energy terrain as illustrated in 
Fig. 14. The heights and depths of these additional barriars 
and mells respectively are listed in Table 2, Fortunately, the 
dynamical investigations reported in this thesis are at such 
high energies that these additional features are not likely to 
have much influence on the dynamics. Also, the shape of the 
inner repulsive mall is knomn to be more important than the 
additional mails. Nevertheless, their existence reveals tmo 
important factors: 

(l) It is not enough to examine the PE contours in the 
interaction region? it is essential to ensure the smooth' 



(no) 
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Flg.lt Same as ffg.2 on unscatod and ynskawed coordi- 
nates. 



(no) 
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Table 2 . Barrier heights and u/ell depths for different LiFH geometries 


36 


to 








38 


behaviour of the PES at large i = 1, 2, 3 values also as 
the trajectories are invariably started at large distances* 

(2) It is necessary to ensure that the fitting of an 

analytic function to the PE values does not introduce any 

spurious bumps and troughsthat luere not present in the original 

data. In this particular case, for example, the original ab 

initio values uiere not reported in the region ijuhere the new 

hills and valleys have been discovered. This is an instance 

72 

where a constrained optimization may be more useful than a 
simple nonlinear least squares fit of the parameters of the 
chosen analytic function to the data. Clearly, a better analy- 
tic fit is required for the LiFH system and to obtain the same 
is not a simple matter. The problem of fitting initio 
surfaces continues to be a ’bottle neck’ in going from electronic 

structure calculations to dynamical studies and posesa challenge 

73 

to the practicing molecular dynamicist. In this particular 

case, the sudden nature of the PES poses a special problem. We 

have already indicated in Chapter 2, our failure to make the 

LEPS function reproduce the major features of the surface. This 

7 0 

experience is shared by others — for example, Schor et al. 

7 4 

in their attempts to fit the collinear BeFH surface. Efforts 
in this laboratory to use a Rotat ed-Morse-Curve-Spline and its 
variations to fit the LiFH surface for the collinear configura- 
tions have not been successful and we attribute it to the sudden 
nature of the PES. 
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Since we are investigating the dynamics in collinear as 

well as noncollinear geometries it is essential to examine the 

angular variation of the LiPH potential in detail, A traditional 

way to do this is to make fixed-bond-distance polar plot in R, fl 

space where R is the cent er-of-mass separation between Li and FH 

and Q is the angle between R and Such a plot is shown in 

Fig, 15 for HF fixed at its equilibrium bond distance and 8 is 

varied. Since we are reporting the effect of reagent vibration on 

reaction dynamics, we are including the polair plots of the PES 

when HF is fixed at its inner and outer turning points of v =0, 

1, and 2 states, in Fig, 16. The existence of the spurious 

bumps and troughs is shown by these plots also. It must be 

pointed out that the PES for the rigid rotor HF-Li shown in 

Fig. 15 differs significantly from the potential reported earlier 

25 

by Lester and Krauss, reproduced in Fig, 17, 


All plots shown in Fig, 11 -16 are displays of different 

75 

cross-sections of the CSCM surface and are deficient in that, 
constraints in the form of fixed LiFH and fixed h-ave been 


introduced. Plots of the PE contours in hyperspherical coordi- 
7 6 

nates may be more revealing .but they remain to be plotted forthis 


system. Until such a plot is available, we have to syhthesize 
in our mind the composite four dimensional PES in the (R^, R 2 i 
Rj) space. 



40 



Fig, 1 5 P E contours for tho CSGM surface in polar coerdinates for 


Li + FH (v) , CSCM 



Ftg.tS Sam« as flg.t5 for different UF. 
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3.2 Computational Procedure 

7 8 77 

The QCT technique has been described elseuihere. * * 

The coordinate system and the procedure for initial state 

selection and final state analysis are the same as in ref. 77. 

7 6 

We hav/e used the REACTS program of 3.L. Schreiber ujith appro- 

79 

priate modificationsi The IMSL subroutine GGUB has been 
utilised for the generation of pseudo-random numbers for 

selecting the initial variables. It has been shoujn to have no 

80 81 
nonrandom behaviour. Stratified sampling technique mas used 

for selecting the impact parameters (b) at 0 .5 % intervals and 

on the mhole 300-700 trajectories mere run for each set of fixed 

values of U, R and T, The trajectories mere started mith the 

reagents set apart at 12 %. All other variables mere selected 

7 8 

randomly by the standard procedure. * Hamilton’s equations 

mere solved numerically by the 5th order Adams-Bashford 

82 

predictor and 6th order Adams-Moulton corrector. The accuracy 
of the integration mas checked by energy and angular momentum 
conservation and step-size reduction criteria, 

3.3 Reactive Collisions! Results and Discussion 

3.3.1 Effect of U, T and R on reactipn cross section 

Since me are investigating the LiFH dynamics as a proto- 
type alkali-hydrogen halide exchange reaction and no experimen- 
tal result other than that. of product scattering for the ground 
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vibrational state of the reagents at selected T are available 
for this reaction, vue have chosen to investigate the effect of 
U, T and R on for this reaction in a systematic manner and 
unravel the major features of the dynamics rather than to report 
rate constants etc. 

■ We have computed the 3D trajectories at T =* 8.7 and 
20.46 kcal mol‘^ for v = 0, 3 = 0. For fixed T = 8.7 kcal 
and 3 = 0, lue have increased v to 1 and 2. At the same T, me. 
have investigated the variation of luith 3 for v = 0 and 2. 
The results are displayed in Fig. 18 along mith the Monte Carlo 
error estimates of 2 cr (95^ confidence) value. 


Under the conditions of our study there is no vibration ^ 
threshold fnr the 3D collisions . N evertheless , t her e i s a la^ 
vibrational enhancement of the reaction. At fixed 3- 0, as 
\l increases from 5.87 kcal rnor^ to 17.24 and 28.11 kcal mol \ 
Increases from 0.75 + 0.24 to B.67 * 0.94 and 12.41 .1. 04 9^ 
respectively. In comparison, for v = 0. 0=0 increase in T fra» 

8.7 to 20.46 kcal mol results in an increase in from 0,75^ 

0.24 to 1.76 + 0.40 only. In particular, repartitioning of 
energy from T to V by ^ 12 kcal mol'^ increases by a factor 
of eight. These results can be attributed^^ toJ ( 1 ) the saddie^^^^^^^^^^^^ 
point lying in the exit channel and (2) the sud^ nature of the 
PES as has been discussed in Chapter 2. The lack of a substan- 
tial enhancement of by T can be explained by the incoming trajec- 
tories reaching the inner repulsive ujall. as if in a 'cul de sac'- 
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similar to the light scattering from a mirror kept at<v 90° incident 
angle. Increasing \l means increasing the momentum along the 
retreat coordinate and therefore, helps in crossing the barrier 
located in the exit channel. Also, the trajectory is able to 
•cut the corner' and reach the product channel readily. 

The fact that the origin of the vibrational enhancement 
is different from that of trnslation is evident from the plots 
of the reaction probability (P^) and partial cross section as 
a function of b in Fig. 19 and 20. Increase in uiith increase 
in T arises from increase in. at the moderate impact parameter 
strata, whereas the increase in S^, with increase m \l arises from 
two factors? (D increase in P^ in each impact parameter stratum 
and (2) increase in the b^^^ upto which the reaction occurs. 

The latter would have its implications on the reactio^ attri- 
butes which depend on the orbital angular momentum (ILI = 

^ '^rel 

The effect of R on deserves special attention. First 
of all, very limited amount of data is available on this topic 
both experimentally and theoretically.^^ For the ground vibra- 
tional stat^ of HF, increasing R increases slightly. Even then 
R leads to an enhancement slight larger than T. AS^^/aR 0.2 
(kcal mor') when compared to As/aT = 0.086 %V(kcal mol-' ) . 

For V =2, the dependence of on 3 is entirely different. 
Initially there is a sharp decline in 3^(3); then there is a 
dramatic increase followed by a possible levelling off of 




fit. 


functioo of i«P»ct for w » . 


fiifr«r«nt Ti 
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u/ith increase in 3 upto 15 ujith a minimum occurring around 3 = 5, 
In the 3 range 5-9, reagent rotation is nearly four times more 
efficient than vibration in causing the reaction* 

R ** R 

The plot of P and 2n b.P u_s b in Fig. 21 illustrates 
houj the decline in 3^(3) arises from a decline in P (3) at all 
the impact parameter strata. Although we hav/e the results at 
other 3 values also, we have displayed the results only for 
3 = 0, 5 and 7 for the sake of clarity. The increase in for 
3=7, over and above that of 3 = 0 arises mai nlv from the large b 
collisions and this would be worth investigating further. 

83 

In a preliminary communication of our results (Appendix- 
II) we had attributed the initial decline in to the disruption 
of the preferred orientation as this reaction has the lowest 
when Li approaches FH at an LiFH = 74 , The increase in 3^(3) 
was considered to arise from vibration- rotation coupling result- 
ing in the HF bond stretch and the molecule behaving as if it 

has been vibrationally excited. Such an explanation was proposed 

20 

first by Polanyi and coworkers. We reiterated the same expla- 
nation on the basis of the fact that for v = 0 there is only a 
slight increase in with increase in 3j but for v = 2, 
change is substantial. We realise that this evidence is not ^ 

strong enough to justify this explanation. Nevertheless it 

19 

remains a fact that there is a cooperative effect of v and 3 
on S^. 




far *2 
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3ust as the molecule rotates aiuay from the preferred 

orientation for small 3, it should also return to the preferred 

orientation for a reasonably larger 3 and this is perhaps tuhat 

is responsible for the increase in at the higher 3 values. 

Rotational velocity greater than a particular value has no 

influence on the orientation alignment and the molecule would 

appear as a blur ♦ This would explain the possible levelling off 

of at 3 = 15. 
r 

The orientation alignment is best illustrated through 
‘daisy* plots which give the position of the .three particles in 
the cent er-of- mass coordinates , for coplanar trajectories. There 
are reasons to believe (see later in the chapter) that the LiFH 
dynamics is dominated by planar collisions. Therefore planar 
trajectories would serve as a diagnostic tool in understanding 
the effect of 3 on for this reaction. A representative 
trajectory that is reactive for 3=0, becomes nonreactive for 
3 = 5 and becomes reactive again for 3 = 9 at b = 1 .5 is 
shown in Figs. 22 -24, At 3 = 0, as Li approaches FH, the 
anisotropy of the potential steers it towards the most preferred 
orientation and the reaction occurs as shown in Fig. 22, We 
should add that in these plots we have displayed only the inter- 
action part of the trajectory and omitted the uninteresting part 
of the approach and the departure. At 3 = 5, rotation of the 
molecule leads to an unfavourable orientation and hence the 
t ra At 3=9, by the time the 
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reactants approach each other the molecule has rotated more than 
once and there is a favourable encounter leading to reaction. 

We have made an explicit study of the effect of the initial 
orientation on this reaction and the results are reported and 
discussed later in the chapter. 

The effect of 3 on S for v = 0 is different from that of 

r 

V = 2. One possibility is that the magnitude of is so small 

for V = 0 that the variation in S^(3), except for the overall 

increase in from 3, = 0 to 9, becomes unnoticeable . Another 

possibility is that there is a steady increase in due to the 

fact that there is an increase in the number of product channels 

becoming available due to increase in E. That is the rotational 

enhancement is ‘statistical* in origin. If this is the case, 

the QCT results mould agree mith the phase space theoretic (PST) 

84 

predictions and also the trajectories mould be long-lived, 

22 

As a matter of fact, molecular beam results suggest the complex 
formation. We also find some of the trajectories to be long- 
lived and to shorn multiple collisions. A particular trajectory is 
illustrated in Fig, 25, The declining r|_j;p shoms that the 
incoming Li atom approaches the F end of the HF molecule and 
normally this mould have resulted in the formation of LiF and 
r^P mould have increased suddenly. But in this particular 
trajectory. After reaching a minimum, rj_j^p increases again, 
reaches a maximum and then decreases before it Cets in to o scillat e 
indicating the LiF formation. This is clearly a case of a 
secondary encounter. 
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The statistical behaviour is relfected in the product 

energy and angular distribution also as will be shown later in 

the text* It must be stressed that not all trajectories are 

long-lived. There is only a qualitative agreement between the 

85 

QCT and PST results for this reaction. Therefore uie could 
only say that the dynamics for the HF (v = 0) shows more statis- 
tical behaviour than the HF (v = 2) and that the decline in 
Sj,(T) observed for v = 2 is distinctly dynamical in origin. We 

should also add that studies of the 3 - dependence of for 
1 9 

other systems show that orbital angular momentum also plays 
an important role and that any attempt to predict the 3- depen- 
dence of Sj, should take into account the angular momentum 
factors also. 

3.3 .2 Product energy distribution 

While a study of the effect of U, R and T on 3^, reveals 
the nature of the PES in the entry channel, a study of the 
product energy distribution provides clues to the nature of 
energy release as the products separate. Its dependence on 
U, R and T provides additional clues on the details of the 
dynamics. > 

Product vibrational energy distribution (pud) for v = 0, 

-1 

3 = 0 at T = 8.7 kcal mol , shown in Fig. 26 is nearly statis- 
tical and it remains so, for other values of 3 ( =3, 5, 9) 
reiterating our earlier observation that the LiFH dynamics 
exhibits certain amount of statistical behaviour for HF in its 
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ground vibrational state. Increase in T to 20.46 kcal inol 

broadens the PVD and populates LiF upto v' = 6 (in contrast to 

v* = 2 for u = O) . Increase in v from 0 to 1 under identical 

•*1 

conditions (3 =0, T =8.7 kcal mol ) also broadens the' PVD and 
v=2 leads to a definitely nonequilibrium distribution as shou»n 
in Fig. 27, Increase in 3 for v = 2 further broadens the PUD 

resulting in a distinct bimo dal distribution for 3=9 and 
broadens still further for 3 = 15, It must be added that the 
decline in S^(3) is associated with an increasing population at 
higher v' and the increase in 5^(3) with an increase in the 
population at the lower v* . 

Product rotational energy distribution (PRD) also spreads 
out with increase in v (from 0 through 2) and T as shown in 
Fig, 28, The characteristic feature of this plot is that with 
increase in V, the maximum value of 3'(3^ that is populated 
also increases. An increase in T shiftsthe peak of the distri“ 
bution from 3‘'V 20 to 3'ru30 in addition to shifting 3^_^ from 
49 to 69, This must have its origin in the increased E as well 
as I Lj , For example, the maximum b at which reaction occurs 
increases from 1,28 % to 1.58 % as T increases from 8,7 to 20,46 
kcal mol , Since the PRD is more spread out than the PVD, it 
becomes difficult to comment on the effect of 3 on PRD. 

The simplest way to characterise a distribution is to 
report the first moment of the dependent variable. For v = 0, 

3=0 at T = 8,7 kcal mole \ <^\/*^= 1.09 _+ 0.34, 2.56 + 0,76, 






ri». ta* Product- rotational atata diatcibution for 3 »0, for 

V «fii^ T 
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^ = light) category and have been predicted to shovu the 
L 3 conversion as the initial reduced mass of the reactants 
IS very much larger than that of the products* When L is 
small (as it would be when the product reduced mass is small) 

— ^ ^ —V — ■ 

at D'VO the angular momentum conservation L + 3 = L + 3 
requires that Lo=i.3'', In the case of the reaction (B4), Li is 
not heavy enough to lead to a 10 dZ conversion. Particularly 
for V = 0, the reaction takes place at small b and therefore, 
there is very little correlation between L and 3 as shown in 
Fig. 29. For v = 1 and 2, however, the reaction is dominated 
by collisions at larger b and hence L - 3 correlation becomes 
substantial as illustrated for v = 2 in Fig. 30. An implication 
of this correlation is that 6 = 0° or 180° and this has been 
found to be the case, this condition also implies coplanarity of 
the collision (see later). With very little change in reactive 
b, increasing 3 from 0 through 15 would mean a lesser correla- 
tion between if and 3^ and this is reflected in the correlation 
coefficient changing from 0.57 to 0.48. 

The fate of additional reagent energy in a reaction has 

been the subject of a series of investigations and conversion 

rules like AU— ^ AU'" and AT AT'^ + AR^ have emerged in recent 
88 

years for several bimolecular exchange reactions. We report 
in Fig. 31 the variation of ^ y, and ^T y with V at T = 

A 

8.7 kcal mol” for 3 = 0. At low U, AV ends up predominantly 
as Ar' , But at a higher V, AW — » A\/'' conversion is more 
important. Detailed studies on model systems have shown the 
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AU A\y^ conversion to be less efficient on a su dden surface 
than on a gradual surface and AT to become mostly AT^ + AR on 
the former. We observe the same trend in our system when T is 
increased from 8,7 to 20,46 kcal mol AU /AT = 0.35, AR /AT = 
0.26 and ATVaT = 0.38. The variation of <\/^>, <(r'> and 
with. R is erratic and no clear cut picture of the R dependence 
could emerge. However, for the particular case of 0 changing 
from 5 to 7 (R = 1.81 kcal mol”^ 3.38 kcal mol <(r'^ changes 
from 10.21 to 13.81 kcal mol"^ implying that AR — > AR' j in 
addition, part of \1 and T have been converted to R . 


3.3.3 Product angular distribution 

The product angular distribution (PAD) in Fig. 32 shows peaks 
at an atomic scattering angle of 0° and 180° for v=0, 3s»0, 

T =8.7 kcal mol” . For reaction proceeding via a long-lived complex, 

the PAD is determined solely by the conservation of the total 
angular momentum. This conservation implies for cases in which 
3 a symmetric distribution proportional to l/sin 0 will 


occur. The distribution we have reported is not completely 
symmetrical. The product molecule is preferentially forward 

scatter.ed. The dynamics is not completely statistical because 

^ -1 - 'O' ■ ^ 

T is larger than the well depth (cw6 kcal mol , for LiFH — 
114°) in the PES. At a still higher T of 20,46 kcal mol , the 
peak at = 0° becomes prominent implying that the product 




( POJSJS/^^) C0p/<»SP 
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molecule is more backujard scattered. For u - 2, the peak at 
/q ■>= 180° essentially disappears. That is >the product mole- 
cule is predofninantly backuiard scattered# This is surprising 
as the vibrational enhancenient arises partly from an increased 
b and under these conditions one expects more foruiard 


max 


-.1 

scattering# For v = 2, at T = 3«7 kcal mol $ changes in 
<( 631 ^ are uiithin the statistical error estimates for 0 = 0-9 » 
but a further increase in 3 from 9 to 15 results in a change in 


/0 . > from 92+4,6^ to 107+7.2°. That is, a large rotational 
\ at 

excitation results in' a larger foruiard scattering of LiF. We 
have no explanation for this observation at this stage# 


3 # 3 #4 Planarity of the collision 

When the rotation of the reactant molecule and the orbital 
motion of the incoming atom uiith respect to the rotor are 
coplanar, there is no Gomponent of motion to move the three 
atoms out of the initial plane and the reaction remains coplanar 
An observable of planarity^° is the polarization angle r, that 
is the angle betiueen the initial relative velocity vector and 
. Results in Fig. 33 shouj the distribution of 7 to be 
mostly around 90° ujhich corresponds to planar collisions. This 
is because the centre-of-mass of the rotating HF is nearly at 
the F-end and therefore the orbital motion of the attacking Li 
atom is nearly the motion of Li about F, in a plane defined by 
the relative velocity vecotr and the centre of mass of FH« 




Fig. 33. of i fos airf«r«*t o T for 3**0 
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Since the H atom departs uiith negligible momentum, Li and F 
continue to rotate in the same plane after formation of the neiu 
bond. The shape of the distribution in 7 can be seen to be 
correlated with the L •— 3 ^ correlation or lack thereof. 

3.3.5 Effect of initial orientation 

2 91 92 

In recent years it has become possible ' ^ to select 

the relative orientation of reactants and study its effect on 

various reaction attributes. Although no experimental results 

2 

are available for the reaction (R4), Zare and couiorkers have 
reported recently on the effect of initial reagent orintation 
on PUD for the strontium analog. We have computed at three 
different orientations a (defined in Fig. 34) = 0 ,90 and 
1 8QP for V = 2, 3 = 0 at T = 8^7 k cal mol . In these calcula- 
tions, a u/as chosen at a fixed value UJhile all other variables 
uiere selected randomly from appropriate distri bution functionSf 
as usual. The results are reported in Table 3 and they shouj 
that is the largest for a = 0°, less for 90° and the least 
for 180°, indicating that the approach of Li towards the F end 
of HF is much more preferred over the H end. The product energy 
disposal also shotos a definite dependencB on the initial orien- 
tation of the reagents. As a increases, the average fraction 
□ f vibrational energy, <(fy''> increases at the expense of • 

Effect of changing a on PWD is shown in Fig. 35. As cc changes 
from 0° to 90°, more of higher v* levels become populated and 
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the. PUD becomes broader. At «= 180°, the P\iD becomes 
narroujer, ujith ths ^ level being preferentially populated. 

It is clear from Fig, 36 that the PRD becomes more 
specific for the broadside attack. Though the maximum of PRD 
fora = 0° and 90° remains unchanged at 0 50, for 180°, PRD 
becomes narrouier and the histogram peak shifts back to 3 -25, 

As a changes from 0° to 90°, the peak of the PAD histogram 
shifts from 6^^ 0° to 6^^/^ 80° as shoujn in Fig. 37 . For 

a = 180°, the PAD becomes broader and the product is scattered 
in all the directions. 

We have also investigated the effect of orientation at 
different 3 values, by calculating 50 trajectories in the range 
b = 1-1.5 % at uihich 3 is found to have a large effect, and 
the results are given in Table The orientation that is 

reactive at 3 = 0 becomes nonreacti ve at 3 = 5 and vice versa. 

At 3 - 9, the reagent rotation is found to have very little 
"effect. This is because, half of the rotational period for FH 
at 3 = 5 is about the same as that of the time taken by Li to 
approach FH , Thus, even though Li u/as at the H end initially 
(0=180°), as it approaches FH, it encounters the F end and the 
reaction takes place readily. A similar argument- holds for 
the non-reactivity of thea= 0° orientation at 3 * 5. At 3 = 9, 
the molecule rotates faster and the orientation effect becomes 


minimal . 
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rig. 37. Preauct «rtgui«r tot «• 



Table Reaction probability as a function of the approach angle CC? 
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3,3,6 Comparison ujith other theoretical studies 

A general feature of the alkali- hydro gen halide reactior|s 
is that they involve the dissociation of a covalent bond and the 
formation of an ionic bond. Historically* such reactions were 

23 

explained on the basis of a *harpoon model*# That is> the 
alkali atom pulls the halogen atom by inducing a strong 
coulombic attractive force between the reactants by transferr- 
ing its valence electron to the hydrogen halide molecule. The 
metal atom in effect uses its valence electron as a ’harpoon * 
to pull the halogen atom. The critical separation r^, at which 
this transfer can take place is given by 

2 

IP(M) - EA(HX) 

""c ' 

where IP(M) is the ionization potential of the alkali atom^i 
EA(HX) is the electron affinity of the hydrogen halide and 
e is the electronic charge. Then 5^, is related to by 


Thei ionization potential of Li is known (124,3 kcal mol ). 
The! EA of HF at each v is calculated following the procedure 
usJd by Polanyi et al.^^ The value of r^ is found to be 
1.8^, 2.00 and 2.16 “ft for v = 0, 1 and 2, respectively. The 
corresponding S values are 11.0, 12,6 and 14.7 h compared 
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to our QCT results of 0,75, 8,69 and 12,41 ^ , respectively* 
Considering the simplicity of the harpoon model, the discrepancy 
is not surprising. The reasonable agreement betiueen the tuio 
sets of results at v = 2 suggests that the harpoon mechanism 
may be realistic at the higher v values for (R4), 

The only other QCT result available for the reaction 

(R4) is due to Zeiri et al,^^ on their semi- empirical ZS 

surface. They have computed an of 2.46 at T =23.3 kcal 

mol”*'. This is larger than ujhat one mould expect from the 

experimental value of = 0.94 % at T = 8,7 kcal mol and 

using the experimentally determined AS^/AT = 0,025 « /(|<cal mol ) 

in the T range 3-8.7 kcal mol’*'. It is also larger than mhat 

, : >1 

uie mould predict based on our QCT result at T — 20,46 kcal mol 
and our computed AS^^/AT in this T range (see above) . This discre- 
pancy can not be due to the difference in for the tmo diffe— 

■ - 1 

rent PCS, The ZS surface has an mhich is '^4 kcal mol 
higher than that of the CSCM surface and therefore one mould 
expect a lomer value of on the former than on the latter. 

The pronounced backmard molecular scattering predicted by Zeiri 
et al.^^ is contrary to the experimental result. Therefore, 
there must be some other topological difference betmeen the 
tmo surfaces. Visual examination of the PE contours shoms that, 
the ZS surface is not as sudden as the CSCM. Collinear QM and 
QCT calculations^** do not reveal any vibrational threshold 
reiterating our viem that their ZS surface is lacking in its 
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sudden character. This would explain the larger value of 
at a higher T on the ZS surface. 

3 .3 .7 Comparison with experimental results 

In this section we compare our 3D QCT results with the 

experimental results available for the LiFH and related systems. 

2 

Our calculated value of = 0.75 + 0.24 % compares well with 

22 Q 2 

the only molecular beam result of 0.94 h at T = 8.7 kcal 

- 1 

mol . The experimental PAD shows a prominent forward mole- 
pular (LiF) scattering in agreement with the QCT predictions 
in Fig, 32, Our calculated <(t''/^ is”^ 60Z of the total energy, 
in excellent accord with the experimental value^^ of "Ct^^ 55^ , 

The calculated distribution of shows prominent peaks at 

0 ^ 1 ^ = 0 ° or 180° with a minimum around 90°, This type of dis- 
tribution is characteristic of reactions proceeding via coplanar 
configurations, and the experimental results also suggest 
that the reaction (R4) is mainly coplanar. 

The excellent agreement between the QCT and the experi- 
mental results noted above implies that the CSCM surface is of 
reasonable validity. Although no other experimental result is 
available on this particular system, we can compare our predic- 
tions with results on similar systems. , The marked vibrational 
enhancement predicted by us is in keeping with the vibrational 
enhancement observed * * for reactions (R 8 )» 

(R21) and (R22): 



Ba + FH 


BaF + H 


(R21) 


Sr + FH > SrF + H (R22) 

In particular, \J has been shoujn to be an order of magnitude 

more efficient than T in causing the reactions (R7) and (RB) , 

The marked vibrational enhancement observed for (R9) by Bartoszek 
17 

et al . suggests the PES of these systems to be of the sudden 
category lending further credence to the quality of the ab initio 
^surface for (R4). 

The decline in 3^(3) folloujed by an increase and a possible- 

levelling off of 3^ vuith increase in 3 has been found expari- ‘ 

19 

mentally for a variety of systems. In particular, Blackiuell 

20 2 "I 

et al. have found such a behaviour for (R9) . Disport et al. 

16 

and Heismann and Loesch have observed a steep decline in 3^(3) 

93 

for (R7) and Man and Estler for (R22). The effect of 3 on 
Sj, for (R4) remains to be verified. Here, uie uiish to stress 
the contrasting 3-dependence predicted for v = 0 and u = 2, 

The effect of relative orientation of the reactants on 

2 

PVD has been studied by Karny et al. for the reaction (R22) 
and their results are in qualitative accord with ours for (R4) . 
That is, a broad side attack populates the higher vibrational 
states of the products more. 

On the whole, one could see that based on the satisfactory 
agreement between the QCT results on the CSCM surface and the 
only available molecular beam results we have ventured to predict 
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a variety of o&servables for the prototype alkali- hydrogen 
halide exchange reaction and uje do hope that this stimulates 
further experimental study of this reaction in the near future. 

3*4 Nonreactive Collisions 

We have investigated the effect of v, 3 and T on the 

vibrational energy transfer (VET) processes occurring in this 

potentially reactive system. The purpose is tujo-foldJ one is 

17 

that experimental studies like chemiluminescence depletion 

assume that the depletion is essentially due to reaction and 

this remains to be verified. Another is that there is a lot of 
94 

interest on the effect of v, 3 and T on the VET processes. 

The total vibrationally inelastic integral cross sections (Sy) 
computed by the QCT method for Li 4- FH collisions are listed 
in Table 5 for different v, 3 and T, along luith values under 
the same conditions. It is clear that 3^ is an order of magni- 
tude smaller than implying that the error introduced in 
neglecting VET processes in the chemiluminescence depletion 
studies is small, particularly in comparison to the vibrational 
enhancement and the rotational inhibition. The values for 
this system are unusually loui - mainly due to the sudden nature 
of the surface, that is, V Is lueakly coupled to T. For the same 

reason, S is nearly independent of T in agreement with the 

94 

results of elaborate QCT studies of Thompson on HF-Ar. It is 
tempting to state that varieb with 3 in the same fashion as 
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S^. But the statistical error estimates are of the same magni- 
tude as and therefore uie are not able to make any definite 

statement on the effect of 3 on UET processes. This is surpri- 

94 

sing as Thompson has observed dramatic changes in with 
changing 3 for Ar + FH collisions. 


We have also investigated the rotational energy transfer 

(ret) processes in Li + FH (u = 2, 3 = O) collisions at T = 

- 1 

8,7 kcal mol . The purpose of this study is to examine the 
effect of Li as a collision partner ujhen compared to the poten- 
tially nonreactive Ar, for example and also to test the validity 


of the poujer gap lam for RET processes as it has been shomn to 

95 

fail for rigid rotor CO 2 -H 2 collisions (Appendix III), Unfor- 
tunately, the Li-FH collision is rotationally inelastic upto 
^Li- FH ~ ^ ^ therefore our studies, on this topic are not 

complete. We report, homever, '\3^(3^+l))> as a function of b 
in Fig, 38. We hope to complete the calculation of ^ 

for RET processes for Li-HF (v = 2, 3 = O) and also for Li— rigid 

96 94 

HF and compare them mith experimental and theoretical results 


on HF “Ar and related systems in the near future 
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CHAPTER FOUR 

DYNAMICS OF A MODEL SIX-CENTER EXCHANGE REACTION 


We model a six-center exchange reaction as occurring 
betiueen a van der Waals (\iDW) molecule* A2““'”B2 ^ diatomic 

molecule C 2 as shomn. 


? C 

! + t ■ ■ AB + AC + BC C R23 ) 



Since this is the first time a QCT study of this type of 
reaction is being reported, UJe describe here the method of our 
calculation in detail. 

We assume that the motion of the atoms can be described 

97 


by Hamilton’s equations* 
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oPi 

' . M_ (i = 1, 2 3 N) (1) 

ujhere are the Cartesian coordinates and the conjugate 

mementa of an N-particle system. The dot on the v/ariabies 
represents the derivative ujith respect to time. H is the 
classical Hamiltonian of the system which for a conservative 
system is defined by 

H = T + U = constant (2) 

where T is the kinetic energy and U is the PE of the system. 
Now the QCT technique involves setting up of the Hamiltonian 
and following the time evolution of the system by solving the 
Hamilton's equations numerically till it reaches a final state 
as specified by the final conditions. 

It has been found that Cartesian coordinates are the 
rost suitable for defining the. Hamiltonian for the systems 
with more than three atoms. Ttie coordinates used in this 
study are! 

(i) molecule-fixed coordinates xyz for A 2 » x'y'z'for 82 
and x“y”z"for 82 * with the C.M. of the respective molecules as 


the origin, 
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(ii) space-fixed coordinates X*Y*Z', with the C.M. of 
the tetra— atomic species (UDW molecule), as the origin^ 

(iii) space- fixed coordinates XYZ with C.M. of the whole 
system ( hexa- atomic) as the origin. 

Using the Cartesian coordinates described in Fig. 39, 

^ i=1 "’i 

M = \/(q^, q2» — -- Q-^ 3) 

The equations of motion are: 


-aw 

35- • 


where m^ = m2 = ... - 


my = m0 = . . . .= m^2'“ ‘”3 


'”13~ '’^14 ~ • • • “ e '"c 


Wo and 


mn, mo and m^ are the masses of the atoms A, B and C respectively. 
In order to solve these equations we need to know the derivatives 
of the PE with respect to the coordinates, in addition to the 
initial conditions (coordinates and momenta). 
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(PlP 2 P 3 ^ 


(X1Y1Z1) (PyPaPq^ 


(X3Y3Z3) 

(P. Py Pz > 

H ^3 ^3 


Cp 4 P 5 P 6 ^ 


(X2Y2Z2) 

^^X2^Y2*^Z2^ 


(qi0qi1^12> (P^P-P_) 

(Pl0PllPl2^ ^ ^ 


(X5Y5Z5) 




^ Pi 3'^1 4^15^^ '^"^6^17 ^18 


(P 13 P 14 P 15 ) ^Pl6Pl7Pl8^ 


Fig* 39 . Coordinates and their conjugate momenta 
in the spaGe''fixed frame XYZ. 
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4,1 Potential Energy and its Deriuatives 

99 

Folloujing the treatment of London, Eyring and coujorkers 

derived the secular equation for the six-electron problem. 

Thompson and Suzukama^"^ as mell as Dixon and Herschbach have 

used this approach to study the 6-C reactions. Their calculated 

5 0 • 1 

(^ 88) u/as higher than the afa initio value of 69 kcal mol . 
In our present study lue luill use the same formalism, but luith the 
inclusion of a sato parameter S, to scale douin the to that 

of the ab initio value* 


The fifth order secular determinant given by Eyring and 
99 . 

couJorkers is of the form 




( 5 ) 


ujhere H. . are the integrals 
^1 




and S. - are the overlap terms 
^3 


®ij l^j/ 


The individual terms of the determinant have been given in ref. 99. 
These consist of Coulombic integrals and exchange integrals 

3. . . which are functions of the corresponding (one of the 15) 
internuclear distances^ Instead of solving these integrals 
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we have made use of the LEPS formulation to evaluate these 
integrals. Thus 


and 






«l .i* "ij 
1 + s 

Q . . “ 3 . - 
J-.l 

1 - S 


( 6 ) 


( 7 ) 


where ''e- • and ^E- ^ are the singlet and triplet energies for 
the diatomic molecule given by. 


-2 [-SijCrij - 


( 8 ) 


^D. . is the dissociation energy (including zero-point Bnergy)t 
r® j is the equilibrium internuclear distance and is the 

Morse parameter. The triplet state curve for the atom pair ij 
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IS repr 


esentsd by the Pedersen-Porter functions: 


‘e;_ = ^D. j I'bxp [-2Plj<rij-rij)] 2 «P 

( 9 ) 


for r* - ^ R and 


^E. - = C. .(r. . + A. .) exp (-O- . r. .) 

iJ ^3 ij iJ 


( 10 ) 


r- • ^ R_* 
1 j c 


for 
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We make use of the parameters computed by Raff et al. 
for H 2 » as these values had been used in a number of trajectory 
calculations for three** and four- atom systems* for mhich the 
agreement ujith the experimental results ujas satisfactory. The 
sato parameter S, mas varied till the for the regular hexa- 
gonal geometry calculated using this method agreed u/ith that of 
the Jib initio value. These parameters are given in Table 6. 

Once the matrix elements are evaluated, the secular deter- 
minant is solved for the lomest eigenvalue using the NICER sub- 
routine package^ mhich makes use of the Householder method 
for diagonalization. The PE for the regular hexagonal 

obtained from these calculations are shown in Fig. 40. The 

-1 

for the hexagonal configuration is found to be 69 k cal mol « 

It should be noted that this calculation predicts that 

2 H 2 + > H 2 + 2 HHV (R24} 

is a possible mechanism for the H 2 exchange reactions, since 
this energy is less than that of the dissociation energy of H 2 * 

An important and the most time consuming process in the 
QCT technique is the calculation of the PE derivatives with 
respect to the various coordinates. This is a straightforward 
process when the PES is known in the form of an analytical 
function. In the present case, however, PE is calculated by 
solving the fifth order secular determinant for each nuclear 
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fig. *0. rtyultt h*K«gon«l »« » function of ^rtri* 
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T able 6 . Potential- Energy parameters 


^D(kcal mol 

r®(A) 

^D(l<cal mol 

Bi'i-') 

C(l<cal mol ^ 

A(%) 


109.46086 
0.74127 
1 .97357 
45.35618 
1.88999 
1113.59090 
0.529167 
3 .166560 


0 . 84667 
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configuration as (see above) the louiest of the five eigenvalues* 
This exercise has to be repeated for all the nuclear configu- 
rations during the course of a collision* We calculate the PE 

103 

derivatives by making use of Hellmann- Feynman theorem (HFT). 


by 


According to HFT the force acting on a particle is given 



( 11 ) 


ujhere'^is the wave function of the system and depends on the 
space and spin coordinates of the electrons* denoted by r* and 
on the N parameters which in turn depend on the nuclear 

configuration. X is any parameter on which H is dependent. 
However, HFT assumes a knowledge of the exact wave function, 
Hurly^*^^ has given the conditions under which the theorem holds 
for an approximate wavefunction also. If F(X(j, Yq., 
is the electronic energy for the particular nuclear configura- 
tion denoted by (Xa, Y^, Z^ ) force = 1, — - m (m nuclei), the 
force (PE gradient) in the x direction on nucleus a is given as 


t r ^ (12) 

. ' * SXa aXa 

The term ^e/^Xce takes into account the direct dependence of E 
on X(j , i.e., the occurrence of X^^ in the Hamiltonian. The 
second term takes into ac count the dependence of on Xq through 
the parameters lovolcing HFT means calculating only the 
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first term. Thus HFT holds mhen the second term vanishes. This 

requires for each 'r, either be independent of X^j t or 

^E/^C =0* The latter holds when ■¥ C ? is chosen uariationally • 

n n j 

In general, HFT holds for an approximate function in ivhich all 
parameters ujhich change with nuclear displacement have been 
chosen vari ationally . 

105 

A generalisation of the HFT is as follows: the secular 

equation can be written in the matrix form as 


H C = E S C 


(13) 


IfM = H-ES, then 


M C = 0 


and also C M = 0 


(14) 

(15) 


Differentiating Eqn. (14) with respect to X, we have 




(16) 


Multiplying by C"*" and using M = 0, we have 


^ = 0 


(17) 


But 


15R _ Sh t ^ 1 
M “ ax “ ax “ ax =* 


( 18 ) 
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Hence Eqn, ( 1? ) becomes 




^ H -A 

^ ax ^ 


E c'*' C 


3 Q 


= 0 


or 


■&X 


~+ hH 7^ T- n+ T. 
L» ^ ^ L •" L. Lr "^rtr ^ 

t* sc 




(19) 


( 20 ) 


If X is taken to be the ith internuclear distance and 

zero differential overlap of atomic orbitals is assumed. 


' bl 

a-i 


^ c 
c'^ s c 


or 




® "^k 


( 21 ) 


uihere "C. is the eigenvector corresponding to the eigenvalue 

K 

E. , uihich is, the loujest eigenvalue in our case. Hence by knour 

K 

ing the eigenvector corresponding to an eigenvalue, it is 
possible to calculate the derivative of this particular eigen- 
value. In our program the diagonslisation is carried out using 
NICER subroutine package, mhich produces appropriately 
normalised eigenvectors, that the denominator 

in Eqn. (21 ) becomes unity* RecefAly, uie have learnt that 
Kroger. and Muckerman had used the same HFT approach in calculating 
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the PE derivatives in their QCT study of hot atom reactions. 

^.2 Calculation of Initial Coordinates and [Momenta 

Initially the molecule A2 is aligned along the z-axis with 
its C.M. as the origin. Hence in the molecule- fixed system the 
coordinates of A2 are given by 

X-] == ~ ^^2 ~ ^ 

z^ = r^/2 

^2 “■"^1/2 

where r-j is the internuclear distance of A2. 

To calculate the initial angular momenta, the rotation of 
the molecule is assumed to be in a plane perpendicular to the 
xy-plane and if this plane makes an angle with the xz-plane, 
the components of the angular momentum areS 


(23 ) 



( 22 ) 
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2 

and U(r,0) = U|^^(r) + ' 

2 

and UjyjCr) = D^[ 1 - exp [“^(r-r^)]] . 

and are the inner and outer turning points. 


From Eqn.(25) and (26) u/e have 



mill be negative in one half of the vibrational phase 

and positive in the other half. The initial Xj y» z components 
of the total momentum conjugate to the molecule- fixed coordinates 

are, 





= m^z^Wy (A 2 ) 

= -m-jZ<jO)^ (^ 2 ) 

= Pv 

= '"l^2“y ^^2^ 
= -m^Z2C0^ (A 2 ) 



(28) 


( 29 ) 
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Having determined the coordinates and momenta for a molecule 
aligned uiith the z-axis, me can nom calculate the initial 
coordinates and momenta of the molecule in any position in 30 - 
space by multiplying by the appropriate combination of rotation 
matrices. The tuio rotation matrices used for this purpose are 



cos S 

sin 0 

o1 

»■ 

Rj(e) = 

-sin 0 

cos 0 

0 . 



0 

0 

1 



cos 0 

0 

-sin 0 I 

Ry(P) = 

0 

1 

I 

3 


sin 0 

0 

cos 0 


( 30 ) 


Similar relations hold for 62 and C2 also. 

The molecule- fixed coordinates (x, y, z) of A2 and 
(x', yS z*) of 82 are transformed into space-fixed coordinates 
(X', Y', Z‘ ) by the follouiing transformations (see Fig, 4 l)S 


^i 


p-l sin cos ' 


^i 

y: 


sin S^ sin 

+ R^(|^-,)Ry(8^) 

^i 



p^ cos 84 

' ' :j 

^i 


(i = 1, 2) 

( 31 a) 
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'^i' 


-p2 sin 64 cos 04 




== 

-P2 sin 84 sin ^4 

f R^(? 2 )Ry( 82 ) 

^1 

.^i _ 

'i 

-p2 cos 64 


^i 


(i = 3, 4 ) 

( 31 b) 

0| and 0 ^ are the orientation angles of ^ 2 * ^2 

orientation angles of B2 in the molecule— fixed frame x‘y*z’» 

0^ and 04 are the orientation angles for the \iDW bond, of 
distance p • is the distance betujeen the C.M. of A2B2 ujith 
that of A2 and P2 is the distance betiueen C.M. of A2B2 ujith that 
o f B2 . 

The corresponding momenta in the space-fixed frame luill 
be given by 






= Ry(e^) 


1 _ 


p 

L 


(i = 1 , 2 ) 


... 


— 

^x. 

X 


p ' ' 

I 




Pz. 

1 

^ ^ ^ \ -r-. - . ‘ 



(i = 3, 4) 


( 32 ) 
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To study the collision betujeen C2 it is 

convenient to transform the coordinates X'Y'Z'to XYZ, For 
this* it is assumed that the relative velocity U = v - C, luhere 
V is the C.M. velocity of Cj and C is the C.M. velocity of 
(A2B2), lies along the Z-direction. If the C.M. of (A2B2C2) 
is taken as the origin of the XYZ frame, u^e have 


2 (m-j + m2) C + 2 V = 0 


( 33 ) 


lijhere m2 and m^ are the atomic masses of B and C respectively. 
Hence it follows that 



The position and momenta of the C2 molecule in the XYZ frame 
(Fig. ‘^ 2 ) areS 




06 



mgb m3(Rs-b"^ 

+^2 ♦m3 ’ my+m2 *m3 


X,Y, Zi) 


(X2V2Z2) 


.(X*Y*Z4) 


(Xe Yg Zg) 


(rnTTirij ♦m3) ’ (m^*m2*i^3) 



Coordinates of the iM>J,oOille« ?i^» 63 *”** ^2 
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uihere is the distance betiyeen the C.M. of (A2B2) and C.M. 
of C2 and b is the impact parameter. The position and momenta 
of C2 in molecule- fixed frame (x''y''z") are calculated as in the 
case of ^ 2 * ®3 ^3 orientation angles of C2 in 

x”y"z’‘coordinates. In the above transformation it is assumed 
that the molecule C2 approaches (A2B2) system from the -Z 
direction and the C.M. of (A2B2) and C2 lie in the YZ plane. 

The initial coordinates and momenta in the XYZ coordinates are* 




mj< - 

^ ^ "’2 


ujhen i 
luhen i 


1 and 2, and 
3 and 4, 


Hence by specifying r- , v. , 3., 'H. , 0. (i = 1 to 4) , 

( i = 1 to 4)^ b and U the initial conditions for 

the system is completely defined. 


4.3 Selection of Initial Conditions 

In the present study lue have restricted ourselves to the 
fixed values of though provision for select- 

ing these variables from appropriate distribution functions 
exists in the softujare package developed by us. The orientation 
angles have been selected randomly as? 

e = cos””' (1 - 2 £) 

^ = 2 n j£| 

T) = 2 7 T isj (37) 

liihere £ is the random number in the interval [-It 3 generated 

80 

from the IMSL subroutine GGUB. The vibrational phase u/as 

selected as for a rotating Morse oscillator and is described 
81 

elseothere. 


4,4 Solution of the Hamilton* s Equations and the Product 
Analysis 

Once the initial conditions are selected and the corres- 
ponding coordinates and rooBionta oaiculated, the equations of 
motion (4) are solved numerically ttsing the same predictor- 
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corrector method used in the REACTS program mentioned in 
Chapter 3. We only mention that me have made use of a step size 
= At = (0.5 -1,0) xIO s. The accuracy of the integration 
mas checked by the energy and angular momentum conservation, 
back integration and stepsize reduction criteria. 

In this section me mill discuss the criteria me have used 
to identify the products and analyse their characteristics. 


In the collision betmeen (A 2 B 2 ) and C 2 the possible 
channels areJ 



A 2 + B 2 + C 2 

(R25) 

?> AB + AC + BC. 

(R26.) 

— ^ 2 AB+ C 2 

(R27) 

■ ^ 2 AC+ 82 

( R28) 

— 2 BC+ A 2 

(R29) 


(R25) is a non-reactive process and (R26) is a genuine 6 -C 
reaction. In (R27), (R28) and (R29), C 2 , B 2 and A 2 acts as 
the energising partner in that order. That is, ittransfers its 
energy to the other collision partner mithout undergoing any 
bond breakage. All the possible combinations of the atoms along 


mith their identification labels are shomn in Table? . To 
identify mhich of the processes has occurred me make use of the 
distance criteria. We ctsuipore all the 15 internuclear distances 



110 


T able 7 , Species labels for different reactive channels 


Process 


Non- reactive 


(R25) 


Six-centered reaction 
(R26) 


C 2 acts as an energising 

partner 

(R27) 

82 acts as an energising 

partner 

(R28) 


Combination of atoms in sets 
of three 


12, 

34 , 

56 




AA, 

BB, 

CC 




13, 

25, 

46 

13, 

26, 

45 

AB, 

AC, 

BC 

AB, 

AC, 

BC 

14, 

25, 

36 

14, 

26, 

45 

CD 

AC, 

BC 

AB, 

AC, 

BC 

15, 

23, 

46 

15, 

24, 

36 

AC, 

CD 

BC 

AC, 

AB, 

BC 

16, 

23, 

45 

16, 

24, 

35 

AC, 

AB, 

BC 

AC, 

AB, 

BC 

13, 

24, 

56 

14, 

23, 

56 

CD 

AB, 

CC 

AB, 

AB, 

CC 

15, 

34 , 

26 

16, 

25, 

34 

AC, 

BB, 

AC 

AC, 

AC, 

BB 

12, 

35 , 

46 

12, 

36, 

45 

AA, 

BC, 

BC 

AA, 

BC, 

BC 


A 2 acts as an energising 

partner 

(R29) 




Ill 


and continue the integration until only three of those 15 
distances are less than a specific distance R^, after u/hich 
the interaction betujeen the species concerned is negligible. 

Then the molecules corresponding to the three shortest distances 
are identified and their energy analysed. 


Once a product molecule is identified, its C.M, velocity 
in space- fixed frame is calculated as 




m . + m . 
i J 


^ ^i Y 


(Py)! + (Py) j 


m . + m . 

V# 


(v;j)z 


(Pz')i+(PPj 

m - + m - 
i j 


(3 8) 


Then the components of the linear momentum in molecule fixed 
coordinates are 


P' = 
P ' = 


P ’ 

^Y, 


P ' 
z. 






(k = i,j) 


(39) 


To calculate the intprpal 




■ momentum, the space- fixed co- 


ordinates are transformed iftfcp molecule- fixed coordinates asl 
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\ ■ 


= Y.^ - 


(m. + m. X.) 


yk = ''k 


-k = \ - 



( m • 4- m 

r 


(m. 

Y.+ n,, 

Y 


! 

tmpTj; 

r 


(mi 

Z.l. m. 

z 

il 


j'' 


= i> j) 


( 40 ) 


The final internal angular momentum of the species (i j) 


=■ 


(y^Pj -z.P ' ) + (y .P ' -z .P ' ) 

2i» X y * ' '*'1 z * j ^ 


= (-iP;.->CiP'.) * (zjp;,-xjp',) 


(m.' .) ' 

■ lyz 


(x.P; -y.P ' ) + (x.P' -y .P ’ ) 
1 Vi 1 Xj_ 1 y . J X - 


The final rotational energy is given by 

/ nii ® \ 2 


'^rot = 


(M. .) 

LjjI 

21 


u,herp = (n;j)2 * , 


I , 

fTi • m * 

.. 1 J 


m. -fm , 
1 J 


and r. j 


r ** « I . 
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The internal kinetic energy of the species (ij) is 


= ■? (P'x" ^ * Pz >/'"i 


(P’^^ + + P*^^) /mj] 

11 1 


and the vibrational energy is calculated from 




U. . is the Morse potential for the diatomic species (i j) . Thf 

final translational energy is calculated by substrating 

and E' . from the initial total energy, 
rot 


4,5 Preliminary Results 

We report here the preliminary results of our investigar* 
tion using the program developed by us for studying the 6-C 
processes. The value of UDW bond distance + P 2 ) has been 

t ak en^ to be 3.0 ^ and R as 6.0 As our objective is to 

explore the conditions under uihich the 6-C pathway is dominant, 
we have restricted ourselves to follow the collision at fixed 
impact parameters. Further, because of the enormous amount of 
computer time taken for the completion of individual trajec- 
tories. we could compute only a few test trajectories. Also 


tories, we could compute only a few test trajectories, f 
i ' ' A/’ ■■ . . _ 

despite the integratipn’civer 10,000 steps t1 step size - 


1 X 10*“'° s ) mo si ^^ori 6s at large imp act p ar amet e r s 1 
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do not go to completion. These trajectories take<x#9 min. of 
cpu time each, in DEC- 1090 computer. 

To illustrate the accuracy of the integration and to 
demonstrate the utility of the developed computer program, me 
give the initial and final coordinates and momenta as mell as 
the results of the back integration for a typical trajectory in 
Table 8. Several test runs mere made for a mide range of 
initial conditions: v^ =0-4, 3^^ = 0-4 (i =1-3) and T = 

10-100 kcal mol”^. 

In Table 9, me are presenting the results for 52 traject- 
ories mith all three molecules being in their ground vibrational 

and rotational states and mith a relative translational energy 

-1 * 
of 80 kcal mol . Under these conditions , all the trajectories 

lead to the decomposition of the UDW molecule as expected. As 

a result me could only examine the energy transfer processes 

occurring in these collisions. It can be seen that T ■— is 

the major energy transfer process occurring in this system. In 

all ‘cases, the molecule-3, i.e. the molecule that collides mith 

the VDW molecule, takes amay the major portion of the energy 

that is transferred. 

In spite of the fact that T is greater than no reac- 

tion has been observed. We should stress that, because of the 
small number of trajectories coin|j«ited, no conclusion can be 
dramn at this stage. An undeiss^^i^ng of the dynamics of the 
6-C processes should study, -using our ^ nem 
programi 
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T able 9 , Average energy transferred (in kcal mol ) at different impact 
parameters 
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CHAPTER FIUE 


SUMMARY AND CONCLUSION 

Our investigation of the dynamics of collinear reaction 
(R4) on an ab initio surface reveals a sharp vibrational thres” 
hold equal to E^ and this uie attribute to the sudden nature of 
the PES. Our 3D calculations do not shoui such a, threshold, as 
the surface is less sudden for noncollinear geometries and the 
sideways approach is the (energetically) most favoured configura- 
tion for this reaction. Nevertheless, there is a substantial 
vibrational enhancement in qualitative agreement with the experi- 
mental results for other alkali-hydrogen halide exchange reactions 
and with other theoretical results on PES with the barrier crest 
in the exit channel. 

Our collinear QCT results' on the CSCM surface agree well 
with the corresponding QM results except, for the oscillations 
in the P^ versus E curves 'of .th® Our attempts td 

construct ' an' LEPS surface wit^.i§^"-%aic^teri sties of the 
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ab initio surface have not been successful in that our best LEPS 
surface is not sudden enough and does not shoiu any vibrational 
threshold. 


Product angular distribution predicted by our 3 D QCT study 

- 1 

for V =0, 3 =0, T =8,7 kcal mol is in qualitative agreement 
uiith the molecular beam results lending credence to the validity 
of the ab initio surface. Based on this success, lue have 
ventured to predict other observables for this reaction. We 
have paid special attention to the effect of 3 for this reaction. 

An increase in 3 (from 0 through 9) for v=0 at T=8,7 kcal 
mol results in an increase in in qualitative accord uuith 
the PST prediction. That is, the rotational enhancement results 
partly from an increase in the number of product channels becoc^ 
ing available mith increase in the total energy. Some amount 
of statistical behaviour for this reaction is reflected in some 
trajectories being long lived and product vibrational energy 
distribution being nearly statistical. 


The effect of R on is definitely dynamical in origin 
for v=2 at T=8,7 kcal mol V, There is a marked decline in S^(3) 
(for 3=0-5) folio wed, by a dramatic increase (for 3=5-7) and 
a possible levelling off of S with further increase in 3 (from 
9 to 15), These results are similar to the effect of 3 on 




observed experimentally for ^ (and alkaline 

earth)-hydrogen^h^#|;-;^|^|ft^^^^|ii#'^,e^..,The ^decline is ’ 
attributed -preferred orientation for 
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the reaction. Oust as the molecule rotates avuay from the 
preferred ‘alignment’ , it can bring the molecule back into 
alignment and this mould explain the subsequent increase. Further 
increase in 0 resulting in increased rotational velocity mould 
only make the molecule appear as a blur and this mould explain 
the levelling off of S^(0) . We must stress that in the 3 range 
5-7 for v=2, the rotational enhancement is nearly four times 
larger than the vibrational enhancement mhich in turn is an order 
of magnitude larger than the effect of T on S^,. 

Our ideas on the effect of rotation being linked to the 
relative orientation of the reactants are corroborated by detail- 
ed QCT investigations of the reastion at cc = 0°, 90° and 180° and 
the effect of increasing 0 at these angles on the reaction probar 
bility at a fixed impact parameter. The initial orientation has 
a noticeable effect on the product energy disposal also. We 
find that for the reaction (R4), a sidemays approach of Li to 
FH results in an increase in the population of the product LiF 
molecules in their higher vibrational states. This result is in 
accord mih the only such experimental result available till 
this date, for its strontium analog. This orientation effect 
deserves further attention. 

In spite of the reaction (R4) being the prototype of the 
alkali- hydrogen halide exchange reactions, Li is not heavy 
enough for this reqetion to the | L { 1 3 j behaviour 

that can be expect^' of this family ot 



121 


reactions, from kinematic considerations. This is reflected in 
the lack of a lOO/ correlation betiueen | Lf and |3'^t for this 
reaction. 

An interesting result of our QCT study is that the distri- 
bution of the polarization angle T is centered around 90'^ indica- 
ting the reaction to be coplanar, in agreement ujith the molecular 
beam results. 

Therefore lue see that the available ^ initio surface for 
the reaction (R4), despite its shortcomings, leads to results 
in reasonable accord vuith the available experimental results and 
has made us go a step further to predict the effect of T, V and 
R on various attributes for this reaction. We hope that this 
study mould stimulate further experimental study of this 
reaction. 

We have compared some of our results mith the available 
experimental results for some of the alkaline earth analogs. 

An important difference betmeen the tiuo families of reactions 
is that configurations such as F-M-H are possible in the case, 
of the alkaline earths and they must lead to certain interesting 
dynamical results. The prototype of this family. Be + FH has 
been receiving attention lately and experimental and theore- 
tical studies of this reaction can be expected in the near 
future . 

■' ?■■■ ' i " 

In Chapter 4 »e tleecrilfe^ ithe methodology, adopted in 
building up ©ur softi^^i^ t© study the six-center (6-C) 
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processes. A special feature of our study is the use of Hellmann^ 
Feynman theorem in extracting the derivative of the potential 
energy. Although lue have not observed any reactive trajectories 
for the model 6-C exchange process that lue had taken up for our 
studies, the preliminary results reveal a dominant T — ^ R process 
in the collision between a diatom and a tetra-atomic van der 
Waals molecule. Apart from opening up the venue for a study of 
several 6-C exchange reactions on model and realistic surfaces, 
our softiuare package makes it possible to study the collision- 
induced- dissociation of van der Waals molecules among other 
things as this topic has been in the forefront in the last 
couple of years. * We hope to get definitive results on 

some of these problems in the near future. 
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The London-Eyring-Polanyi-Sato (LEPS) function is defined 

ast 

3 3 1/2 

r^, t ) = r Q -[^ E J 

1=1 1 j=1 

inhere 

Qi ' i + ^EiCr^)], 

= 2 ' [ r^) - E ( r^) 3 , 

^Ei= D^X^(X^-2), 

^E^= [(l-S^)/(l+Sj^)](D^/ 2 )X^(X^+ 2 ), 

and 

= exp [“0j^(rj.- rp] , 

The spectroscopic constants D^, inhere i = 1-3 are 

dissociation energies, Morse parameters and equilibrium bond 
lengths respecti uelyj are the Sato parameters. These 
constants are listed belotoS 

D (eU) 

r“ (8) 

e ( 8 "^) 

S 


LiF 

HF 

LiH 

6.0030 . 

6.1150 

2.5170 

2.9554 

I.TjJP . 

3.0147 

d''.S366l\;:r; 





4: ‘i4,4 ' 

0.5 1 - 


*’5 
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reliction crt>ss kcctuin 6V^) dccrc.jscs mtii.ilK ,incl thcit 
with an increase in J The lorr.irr effcii i** .ittnbtiied to the 
disrupiHW of the favvwed oricntathin for the ccaction The latter 
effect ts explamcd on the h isis of the T II bond 'itreiching due 
to centrifugal disiori ton iit large {J, v). I ndcr the conditions 
employed in this study, at large J, rt^a^^ent rinatian n neariv 4 
times more efjiaeni than rva^em i ihnitKm. which in turn is more 
effective than reagent translation in causing the rcactum 
Although the last 20 years have witnessed an increase m un- 
derstanding of the effect of reagent translation and vibratum^ t>n 
the rates of chemical reactions, the study of the effect of rcagcni 
rotation has been limUcd, and as a result, the understanding of 
the role of reagent rotation in chemical reactions ha^^ remained 
poor. 

The OCT studies^ of the effect of J on have m«tly focused 
their attention on the reaction 

+ H (2) 

and its isotopic analogues. Depending on ih*: PFS employed, the 
effect of increasing J on S,. was varied; ( I ) a dramatic drop in 
S, foilowcu by a leveling off of the same: (2) a slight increase in 
Sy from J =5 0 to I followed by a decrease in 5,: (31a substantial 
initial decrease followed by an increase in S^, Experimentally, 
Kletit and Persky’ showed that the rate of reaction 2 was nearly 
li^tisttive to J in the range J = 0--2. Bernstein’* and co-workers 
^owed that a small increase in the reagent rotational energy (/?) 
resulted in a small increase in the reactivity for the reaction 
K Hh PCs (KFCsI — K¥ + Cs (3) 

and a small decrease for the rcact^cm . ^ 

k + FRb ^ fKFRb]— Kr + Rh " (4) 

In both cases, the effect of reagent rotation on Sy was anabgous 
to that of relative transtation. ,i 

In recent ycarsl there have been some experimental resniis 
available for some alkali atom-hydrogen halide reactions. 
Blackwell et al.’ concluded from chemiluminescence d'^pktion 
ex^riments on the reactions 

Na -F XH {X * F. Cl) — NaX -F H (5) 

that there was an initial decrease followed by an increase in the 
reaction rale with increase in J. The initial decline was also 
ot^rved by l>jspert et for the related reaction 

K-FaH — KCl-f H (€) 


^fcct of Reofeii Rolatioii an Cross Section Imr tile 
RdK^tiofi U -f FH £ 


L NoorBatcha**** and N. Salhyamurthy(Rt#*‘^ 


Department of Chemistry 
tndian Institute of Technology 
Kanpur 208016, India 


Htme of these reactions has been studied so far theoretically to 
ttii^rstaml why there is a decreasing/ tnercasing effect of / on 
the reaction rate. Recently, howevw, for the simi^est alkali 
atom-hydrogen halide reaction (I), a faffly accurate ah inido PES 
has become available^*^* and: has also bcai fitted lo an analytic 
functioii/ Tliercforc we fptuid this to be an ideal sy^em for which 
the effect of reagent rotation reaction cross section could, be 
studied theoretically. 

We have carried out QCT calculations for this reaction on an 


Received September 9, 1981 

in this communication, we repoh on the effect of reagent 
rotation on the cross section for the reaction 

Li-FFH — LiF + H (1) 

for HF in its a « 2 vibrationat state in tlte range of rotatiofial stale 
0 / S 9, at a relative tramlalkmal energy erf* T * 8.7 kcal ^Krf'^ 

based Oft three-dimensional quasi-classical trajectory^ {QCT) 
studies on an ab initio potential-energy surface {FES)>^ The 

(II (a) fit partial fiiintinieiit of the nxittirem^itt for the degree of Doctor 
of Phftoe^y. (h) On leave from the Amcrtcait College. Madurai 625010, 
Bidia. tmder the FIP scheme. (c> This research was snpported tn part by a 
fraiit from the Department of Sctence and Technology, India, (d) The 
ealct^tions reported tn tins paper were earned out on a DEC- 1090 computer 
at the fndum institute of Teclmology, Kanpur. 

{2| Foruer. R- foi’fk I- M. In *Dywte«n0'of Mdbcirfi»'Colistoiis*;paii« ' 
B; Miller, W. H , Ed,: Pleniim Fress: New York^ 1976; Chap^ 1. 

13) Chen, M. hi. SdMefcf, a K |ff / Chm. 4^ 

14) Carter, S-; M«m^. i. H Mat. #1^^. iWt, 41, 567. -v ? 


ab initio surface. Wc have chosen t? « 2 since the chemtiu* 
ntinesccRce depletion experiments^ on Na -F FTf had v iii the rai^e 
I -6. The value of T - 8.7 kcal mol'* employed in this study is 
the same as that cmph^cd in the only molecular-bcam study*^ 
of this reaetkm (f >. The tktails the QCT method arc described 
elsewhere.^ We mentfon only 'bat the impact parameter was 
samfrfed in a stratifted manner aiKi other varial^ of orientation 
angles and viln'ational phase were sctelcd randomly. 

15) M.; Wolfrum, J. Amm, Rer. Pkfs, CHrm. im if. 47. 

(6) Folanyi, 3, C.; Schreiher, 1. L. Farmtay EHsaiss, Ckem, Soc, 1977. 62, 
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Sc iWtor, A. E; Fope, W. M.: Bcmsteiii, 'R, ft. i, Ckem. Fkys, 
W7,m 34^ T^itdec. tc R, B X Ffm IfTi. 6A 3760. 

IF) m A ; Fiihtiiyi X Cc Smn, 1 1 Ckem, Fhys, I97f, JO, 299. 
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53IX • ’ • 

a Cc Yirdky, ft. N Far&^y IMscms, Ckem, Sac, 
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RIkcoi mof'l 

0-0 07 ? ' 1.81 ^ 3. 38 ' 5 43 , ! . , ■ ' 


2q|._' . ti + FH C V = 2, J) LiF +H 
I T = 8 7 kcQl mof’ . ^ 

18H . ■ " .. 



Hgwe 1, Rmctioa crew sccticm as a fiiisctioii of the rotatioBal state 
Jikmu sette) for HF (p »= 2) at T «* S.7 kca! rod'*. Tlie corespcaidiiig 
nMttoial iroeigief (iis kcal rod'*) axe shown in ^ npper scale, llbe error ' 
ten correspott! to a 95% ccmiideiice icvel. 

The results of the C^JT study arc plotted hi the form of Sj m 
a fiiaaiofi of I k Figure 1. With kcreasc ia J frcrai 0 to 5, there 
is t mouotonlc decrease k S/, further increase k J to 7 aiid 9 
faults in a dramatic increase in 

For reactlcm 1 , ah initk stiMiies^ s^w that an angular ai^Hoach 

li to FH at an LiFH angle of 74® has the lowest terrier, 10 
kcal Ind■"^ fer the reaction. Also, detailed anaiyds of the initia! 
omditiros id tte reactive trajectaies shows that there is a p^ened 
cone of reaction at the F end of the HF molecule. Therefeire It 
is understandable that an increase k J from 0 to 3-5 dmups tte 
preferred oriexitation of the reaction, resulting in a drop in the 
magnitude id In the absence of any additional effect, we wouM 
have expected a leveling off of S,. However, k this paxticukr 
case, b^use of the large vibrationai energy present k the 
mdecule, increase k J results k substantial coi^ugal distialkm. 
This means a stretching of the FH bond under attack, and the 
result would be similar to increasing the vibratimial energy. For 
example, when v is changed from 0 to 1 to 2, increases from 
0J ^ 0.2 to S.7 :h 0.9 to 12.4 ± LO at T » 8.7 kcal md~^ for 
/ * 0. This is explained on the basis of the sudden character of 
the it must be added that such an explanation for the 

observed depciKicnce of on J was suggested earikr by Bladcweil 
et at® although no dynamical results were availaMe at that time. 

At large /, even though the effect of reagent rotation on k 
qualitatively similar to that of reagent vibration, there is a 
quantitative difference between the two. In the energy range we 
have studied, ASJAV = 0.5 A^Ckcal mor*) and AS JAR = 2.7 
AV(hcal mo!”*); that is, rotation is nearly 4 times more eTicient 
than vibration in enhancing the reaction cross section. Reagent 
translatkn is the laist effective in this case with a ASJAT value 
ofO.OS6 AV(kcal mor^). 

In summary our findings arc significant for these reasons: (1) 
Ours is the first dynamical study of the efto of reagent rotaticro 
f« 1 prmotype alkali atom-hydrogen halide reaction, on an ab 
kitio PES. (2) Increase In / results in a substantial decrease 
followed by a rapid increase in with a minimum incurring 
wwnd i » 5. The results arc in qualitative accord with the 
experimental rttults on the related alkali atom-hydrogen, halide 
exchange reactions. (3) For the first lime we have shown that 
reagent rotation can te more effective than reagent^ vibration k 
enhancing a reaction. 

mgrnry n». Lh 7439-93-2; FH. 7464-39-3. 


Cil) Piiaiiyf* J. C.; Sathyamiifiliy* N. Ckem\ 1978, JJ, 287. 
I.; Stitiyarounhy, N, / im ptm. 
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The validity of the power gap law for rolaUonal energy traitifer has been tested in the case of rigid-rotor COj -Hj coiii- 
siom. The power law with a parameter y ta rfiown to be adequate only for iinait amounts of energy transferred. Lor 
larger amounts, uw of a different y seems to be warranted suggestiiv a dcflsdency of the power law or operation of two 
mediMisms for rotational enerity transfer. 


I. Introdttcdofl state-tostate IICS,gf the degeneracy of the fUiabtate 

» 2/f + 1 and 7f the final collision energy, has shown 
^ Polanyl «id Wooddl { 1 } proposed that the proba- a prominent sigmoidal behavior rather than a straight 
biUtyoftraittition from an initial rotational state /{to line. 

a final rotational state Jf decreases exponentially with Brunner and co-workers {12-ISl proposed thaj 

increase in absolute energy gap iA%| = |is>f--Eyjl r m ^ a}AEtfi-y . (2) 

between the two states: *'i"*’'i '• 

P cc exrrf-ri AF rt CD Af « “ 2/f + 1 or iV « /Yq >» (2J^ + i)/(2i, 

*Ap( ClAfcjfl), i ) being the maiier of and. 

where C Is a constant for a given collision system. Fit- with p tlw translational density of states and t a con- 

ting of siate-to-state integral inelastic ertm sections slant. They emphasized [13] “the powe law is su- 

(IICS) and rate cwistant data to the exponential gap perior to the exponential gap law for an of the RET 

relationship |2} has been drown [3] to be equivalent data sets ermsidered when lAE'ifi is less than 20% of 

to a linear surprlsa! analyds 14J . Heller {5) tried to the initial kinetic energy”. With more refined data 

provide a theoretical justificati<tn for such an expo- becoming available over a wide range of /j (= 4, 1 6, 

nentiai dependence. Some computed inelastic cross- 38, 66, 74 and 100) for Na 2 (/i)-Xe collisions, they 

section data have been fitted to the exponential gap found recently 116] that althou^ the power law pro 

relationship (see refs. {6- 1 1 ]). Apart from the com- vided a good representation of the variation of 

paction of data, such an approach has rendered acora- with 1 AF|ff for any given it was not adequate for 

parison with tlic experimentally determined rate data fitting over the complete range of/j probably 

easy [7]. A study ofthe role of the collision partner, because y was not allowed to vary with /j. The energy- 

tbe effect of cfiange in initial ccdlision ener^ Fj and ewrected sudden (ECS) scaling law of DePristo and 

rotational state /•, the validity of classical wedianics P 8| proved to be much better. 

to rotational energy transfer (RKT) procisses* etc. : , Ip we analy^ the st^fe-tb-stalc IICS 

could be directed to a single paratneter Crather than for (;X) 2 -H 2 dotlision# {^J and 

. an array of IICS J 10). However, ttiere arc limitations t^l^icttito the power law with a single ex- 

to tlie exponential gap relationship. In many instances If *1114 c^y up to lAA'jfl/Fj values as low as 

(for example, sec refs. j‘),l I ]) a pitU of the surprlsa! a single-parameter power 

lB(oyj-./f/?rF/^^) versus 1 where is the jatf iwd randoiU; The deviating data 

264 0 00<i 26 1 4/8 i /OWX)- 0000/S 02.50 © North-l iolland Publishing Company 
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points fall on a strai^t line witli a dislinciiy differeni 
'J'hl^’ indicating two different mechanisms for 
the rotational energy transfer ot a deficiency of the 
power law. Wc have also Investigated the dependence 
of'yonJ'j andfj. 


2 . Resuita and discussion 

The qualitative way to check the goodness of fit 
of the power law to the statc-to-statc IlCS is to plot 
against In to see if a strai^t 
line with a neptive slope y results. Such a plot for 
C 02 (/j * 0 )-H 2 (/j “ 0) collisions at 7j * 02 eV fa 
shown in fig. 1 for the two different choices of M It 


is clear iliat the results at small are fitted well hy 
a straight line with a slope of up to A£|f/r} "• 0.08, 

beyond which there are large deviations from the 
St rai^it line. Such deviations were noticed by Pritchard 
et ai. |l3j for the N 2 -Ar system, but only above 
! Ai'jj-1 values whkh were 20% of Tj, They ipored 
such deviations. We find them to be too fatp to be 
neglected in the case of CO 2 -H 2 . As a matter of fact, 
the deviations are larp enough and overt wide enou^ 
A£|f to justify a separate straight-line lit with a slope 
of When ail the available 2 <Jf < 26 
are fitted by a s^le power law, x^/p as defined in 
ref. [ 13 J is equal to 1 7.1 ; when only the first nine 
values of 2 <yf < 1 8, are used, x^f*> drops to 
0 J 8 . Jendrek and Alexander {19] also noticed in the 

case of liH-He coilfaions that the use of a ilnidt)' 



parameter power law Is justified for snnli I only. 
Although they have not mM so, it fa clear tlurt their 
values of ln(c 7 j^/j./ffr/'^) at large )A£{f { do wartent 
the use of a different 7 In that region. Thus it Is evi- 
dent that a sin^e iwwer Itw is not adequate. It fa also 
possible that this arises because of two different tneulia- 
nisms for RET processes. For exmnpk, in the case of 
CO 2 -H 2 It was shown In fig. 7 of ref. {9| that “small 
AJ transitions occur predominantly at large impact 
parameters (b) while large AJ traitiitions tend to occur 
at small to moderate values of 

it may be argued that the deviations of the type 
noticed by us may be accounted for by t^ng tbd ' 
ECS relationship to fit the data. Hie difficulty fa titat 
the scaling applies to de-excitation (/f </|) ptn- 

cessesin]: 


’'(0 ■ P) 

where (: : I) is a 3-/ symbol and Aj} an adiabatidry 
factor, in principle, one could relate the drove equa- 
tion to excitetion cross sections by invokii^ the prm- 
dple of ndcipic^dc reversibility { 20 ] : 

( 4 ) 

Unfortunately, this involves values of <r at different 
coiMshm energies. Also, the IlCScaiculated by the 
quasiclassicat traj«:tory method ( 2 i ] do not satlsiy 


Ffa, }. Power taw fHt to ttate-to^ftate inteersi tnetasiic cross 
sections for COi(/| ■ = 0) eolUsiona at 7j * 02 eV 

using ft • and N*/fo- The arrow indicates apptoximaiciy 
the point at which the switch from one power taw to another 
begini to occur. The'error bars represent to values. 
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t'i|. 3. Dependenc* of on T{ fot difteient dioicet of faiiM rottiioniit itato* of COj «id Mf . 11* ttm cttioutw niHOKat 2® 
valuci. 


nearly the same for boifv. Howe* ?r, In wnne 
C 8 se$, the de-excitatlon cross secttontrfo deviate from 
the straight line(s) much more than their excitation 
analofi. The y values reported in this paper were ob- 
tained mostly by fitting the excitation cross sections 
only. In the particuiar case of^j ■ 30, eq. (3) could 
have been directly tested in piedfcllng O 3 o-*if 
/f < 30 provided o£_ 0 (r|) were known. But these 
values are not known. Substituting for the latter by a 
power law expression as was done in ref. (16] would 
not have been of much use as that will only test the 
power law. 

Values of 7 ,^^ and 7 |,|^ for the different 
and at different fj are listed in table I . The 
values are much more dependable than the 7 |,^ values 
as the former were obtained from a larger number of 
0 values. Therefore we focused our attention on 71 ^^ 
in studying the effect of the initfal rotalionsd state and 
the collision energy. Dependence of 7 |o* on Tj is dis- 
played in fig, 3 for different and.^* . It is seen 
that remains virtually unchanged with chan^ In 
/.“* from 0 to 2. On the other hand, Uierc » a notice- 
able change in 7 j^ withyf (*0, 1 6, 30) although 
not in any systematic marmcr. The value of Tb*‘is 
relatively independent of Tf, a conclusion preylpusly 
arrived at hy Pritchard el al. (13) for N 2 --Mf In con- 


tittt^theexf mthtlgapparameterChadbceadiowit 
to decrease mrHtotmtk^y with increase hi { 9 . 10 . 22 {. 
It must be added that, at Imt for CQj-Me, 7 * 1 ^ does 
aeem to depend on f| {23). 


The power gap hw for rotational energy transfbr 
processes proposed by Brunner et id. { 12) fftt the 
lUte-to-state liCS vidues for ri|^-rotot CO 2 -H 2 cof 
lirions better than the exponential gip law. The dope 
of the straight line for a |dot of ver- 

sus i A£|f| however, changes diarpty rtom to 
TtiM- somethnes at iAffft at tow as 8 % of the hilthd 
comston energy , showh^ the htadequacy of the sbigle- 
Ittrameter power law. Alternatively, the change in 7 
it suggestive of two different mechaititmt for RET 
processes in CX) 2 -H 2 . We do not have sutBdeni in- 
formatkm to distinguish between the two poaibiiities. 

Ttie avidlabte RET data for the C02~^2 system 
surest that y^g^ may depend sli|d>tiy on t^ initial 
lotatkmal state (of the primary rotor) and remains 
unaltered vidfh change in T^. We dso find 
be piefi^dde to iV * Vq for the CXij-Hj system. 
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